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A kinetic treatment is developed for collisionless magnetized plasmas occurring in high- 
temperature, low-density astrophysical accretion disks, such as are thought to be present in some 
radiatively-ineflicient accretion flows onto black holes. Quasi-stationary configurations are investi- 
gated, within the framework of a Vlasov-Maxwell description. The plasma is taken to be axisym- 
metric and subject to the action of slowly time-varying gravitational and electromagnetic fields. 
The magnetic field is assumed to be characterized by a family of locally nested but open mag- 
netic surfaces. The slow collisionless dynamics of these plasmas is investigated, yielding a reduced 
gyrokinetic Vlasov equation for the kinetic distribution function. For doing this, an asymptotic 
quasi-stationary solution is first determined, represented by a generalized bi-Maxwellian distribu- 
tion expressed in terms of the relevant adiabatic invariants. The existence of the solution is shown to 
depend on having suitable kinetic constraints and conditions leading to particle trapping phenom- 
ena. With this solution one can treat temperature anisotropy, toroidal and poloidal fiow velocities 
and finite Larmor-radius effects. An asymptotic expansion for the distribution function permits 
analytic evaluation of all of the relevant fluid fields. Basic theoretical features of the solution and 
their astrophysical implications are discussed. As an application, the possibility of describing the 
dynamics of slowly time- varying accretion flows and the self-generation of magnetic field by means 
of a "kinetic dynamo effect" is discussed. Both effects are shown to be related to intrinsically-kinetic 
physical mechanisms. 

PACS numbers: 95.30.Qd, 52.30.Cv, 52.25.Xz, 52.55.Dy, 52.25.Dg, 52.30.Gz 



I. INTRODUCTION 



This paper is part of an investigation concerning the theoretical formulation of kinetic theory for collisionless 
astrophysical plasmas in accretion disks (ADs) around compact objects, and its application to the study of their 
equilibrium properties and dynamical evolution. Note that what is meant here by the word "equilibrium" is in 
general a stationary-flow solution, which can also include a stationary radial accretion velocity. 

In contrast with the majority of previous treatments, which are based on fluid approaches within the context of 
hydrodynamics (HD) or magnetohydrodynamics (MHD) here we adopt a kinetic approach. This provides a 

phase-space treatment allowing us to formulate a consistent description of plasma dynamics. Kinetic theory is essential 
for studying both stationary conflgurations and dynamical evolution of plasmas when kinetic effects are relevant, such 
as ones associated with conservation of particle adiabatic invariants, temperature anisotropy, finite Larmor-radius 
(FLR) effects (as pointed out in Ref. [6] ) and kinetic trapping phenomena. These properties are relevant for magnetized 
plasmas and in p articular for those arising in ADs [S-Q whenever the plasma is regarded as collisionless or weakly 
collisional [gI, 

In the context of astrophysical ADs, there are several examples of collisionless plasmas of this kind, with both strong 
and weak magnetic fields. One is the case of radiatively inefficient accretion flows (RIAFs) [ill, 113 ' geometrically 
thick disks around black holes consisting of two-temperature plasma, with the ion temperature being much higher 
than the electron one, and the timescale of the Coulomb collision frequency being much longer than the inflow time. 
Other interesting applications occur in ADs around neutron stars and white dwarfs: in the inner regions of such disks. 
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where the magnetic field of the central object becomes dominant, ions and electrons can be collisionally decoupled 
and sustain different temperatures. This happens, in particular, if the radiative cooling time-scale of the electrons is 
much shorter than the time-scale for electron-ion collisions. In this way, electrons and ions are thermally de coup led: 
the two charged species acquire unequal temperatures and the accretion flow becomes a two-temperature flow [l3','T3]. 
In our earlier paper 

(Ref.3, 

hereafter referred to as Paper I), we presented preliminary results in this direction, 
concerning formulation of kinetic theory for investigating stationary solutions for collisionless AD plasmas, focusing 
on configurations with locally-closed magnetic flux surfaces. The present paper is intended as a continuation of the 
previous one, with the aim of generalizing the previous solution to arbitrary magnetic field configurations, which are 
no longer restricted to localized spatial domains in the disk. We refer to Fig.l below and the discussion in Section II for 
an explicit comparison of the two configurations. More specifically, in Paper I the treatment concerned collisionless 
magnetized plasmas characterized by locally closed and nested magnetic surfaces. For such configurations, it was 
shown that suitable kinetic distribution functions (KDFs) are permitted, describing both kinetic and gyrokinetic 
(GK) equilibria (see definition in Paper I), which are represented by generalized Maxwellian and bi-Maxwellian 
KDFs. A main feature was the inclusion in the kinetic treatment of both temperature anisotropy and FLR effects. In 
particular, in Paper I and in Ref.[l5'|, it was proved that these equilibria can sustain a stationary kinetic dynamo. As 
a basic consequence, it was found that both toroidal and poloidal equilibrium magnetic fields can be self generated 
for quasi- neutral plasmas, without ongoing instabilities and/or turbulence phenomena. In particular, in closed nested 
field configurations [and hence in the local absence of net accretion flow], the toroidal field component was found 
to be produced by diamagnetic effects driven by the species temperature anisotropies. As a further development, in 
Paper I and in Ref.[l^, it was pointed out that the kinetic treatment allows one to construct exact fluid equilibria 
(identically satisfying the corresponding fluid equations) . Using a perturbative expansion, a well-defined set of kinetic 
closure conditions was determined analytically for the relevant stationary moment equations. 



A. Accretion disks in astrophysics 

Despite more than forty years of observations and theoretical investigations, there is a lot remaining to be understood 
about the physical processes governing the structure and evolution of ADs. They are observed in a wide range of 
astrophysical contexts |l| and consist of plasma orbiting a central object with the velocities of the inward accretion 
flow usually being much smaller than the rotational velocities. In order for the accretion to happen, there needs to be a 
net outward transport of angular momentum and there are several conceivable mechanisms for producing this (see for 
example pTi. [l8|). The most obvious one is fluid viscosity, but this would need to be an "anomalous" viscosity, driven 
by some type of turbulence, rather than a standard viscosity connected with Coulomb collisions (Spitzer viscosity) 
which would be much too small to explain the observed accretion rates under the conditions actually found in accretion 
disks (see [l9| for a review of turbulence mechanisms). However, other collisionless physical mechanisms are possible 
in principle, such as kinetic instabilities, radiation effects and magnetic reconnection. The aim of the present paper 
is to help in preparing the way for a discussion of these. We focus on AD plasmas immersed in slowly time- varying 
magnetic flelds, characterized locally by open nested magnetic surfaces. For these systems, no kinetic treatment has 
been available up to now. The origin of their magnetic fields varies depending on the type of the central object: in 
the case of black holes, the fields are only ones self-generated by currents in the plasma itself via dynamo effects, 
while with neutron stars and white dwarfs there can also be a magnetic field intrinsic to the central object 0] . 
The interplay between magnetic fields and accretion plasmas can affect the overall velocity profile of the disk, as well 
as giving rise to species-dependent velocities and rotational frequencies j^, i5, JO]. Moreover, the magnetic field can 
be a source of anisotropies in the KDF and allow particular symmetries which influence both the single particle and 
collective plasma behavior. The transport of angular momentum, the accretion flow and the possible generation of 
jets [2Cll422| are all strongly dependent on the magnetic field structure and so magnetic fields play an important role 
for AD physics. 



B. Goals and scheme of the presentation 

The purpose of this paper is to formulate a comprehensive kinetic treatment for collisionless axisymmetric AD 
plasmas including both accretion fiows and collisionless dynamo effects. We include general relative orderings between 
the magnitudes of the external and self-generated magnetic fields and allow the magnetic field be non-uniform and 
slowly time- varying while possessing locally nested open magnetic surfaces. 
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Extending the investigation developed in Paper I, we do this by constructing particular quasi-stationary solutions of 
the Vlasov-Maxwcll equations, characterized by generalized bi-Maxwellian phase-space distributions (sec also Paper 
I), which are referred to here as quasi-stationary asymptotic KDFs (QSA-KDFs). As discussed below, the functional 
form of these solutions is physically motivated. We will show that this makes possible the explicit inclusion of both 
temperature anisotropies and parallel velocity perturbations in the QSA-KDFs (see the definition below in Section 
4). This is done, first, by developing an "ad hoc" formulation for GK theory in the presence of a gravitational field, 
making it possible to directly construct the relevant particle guiding- center adiabatic invariants. The QSA-KDFs 
are then expressed in terms of these. Remarkably, this allows also the consistent treatment of trapping phenomena 
due to spatial variations both of the magnetic field and of the total effective potential (gravitational EM trapping). 
Second, the QSA-KDFs are constructed by imposing appropriate kinetic constraints (see Section 4), requiring that 
suitable structure functions (see below) which enter the definition of the QSA-KDFs, depend only on the azimuthal 
canonical momentum and total particle energy. By invoking suitable perturbative expansions, it follows that the 
relevant moments and moment equations can be evaluated analytically. The solution thus obtained can be used 
for investigating the quasi-stationary dynamics of magnetized AD plasmas, including description of quasi-stationary 
accretion flows and "kinetic dynamo effects" allowing for the generation of finite poloidal and toroidal magnetic fields. 
In particular, the kinetic theory predicts the possibility of pure matter inflows as well as the independent coexistence 
of both inflows and outflows. 

The paper is organized as follows. In Section 2 we summarize the basic assumptions and definitions of the theory. 
In Section 3 we formulate the GK theory for magnetized accretion disk plasmas, deriving the relevant integrals of 
motion and guiding-center adiabatic invariants and discussing the particle trapping phenomenon. Section 4 deals with 
the construction of a generalized asymptotic stationary KDF, with the inclusion of parallel velocity perturbations and 
the adoption of suitable kinetic constraints. In Section 5 we give an analytic expansion for the KDF and discuss its 
main features. Section 6 deals with the relationship between kinetic theory and the corresponding fluid treatment, 
which concerns the validity of moment equations and the analytic calculation of fluid fields. Section 7 is dedicated to 
discussing the temporal evolution of the GK equilibria and derivation of the dynamical equation for the GK KDF. 
In Section 8 we investigate the implication of the kinetic solution for the Ampere equation and the existence of the 
kinetic dynamo effect. Then, in Section 9 we discuss the treatment of quasi-stationary accretion flow within the 
present formulation, showing that solutions with net radial accretion are admitted consistently with the constraints 
imposed by the Maxwell equations. Finally, Section 10 contains a summary of the main results with closing remarks. 



II. BASIC ASSUMPTIONS AND DEFINITIONS 



Ignoring possible weakly-dissipative effects (Coulomb collisions and turbulence), we shall assume that the KDF and 
the EM flelds associated with the plasma obey the system of Vlasov-Maxwell equations, with Maxwell's equations 
being considered in the quasi-static approximation. For definitencss, we shall consider here a plasma consisting of at 
least two species of charged particles: one species of ions (i) and one of electrons (e). 

Following the treatment presented in Paper I, we shall take the AD plasma to be: a) non-relativistic, in the sense 
that it has non-relativistic species flow velocities, that the gravitational field can be treated within the classical 
Newtonian theory, and that the non-relativistic Vlasov kinetic equation is used as the dynamical equation for the 
KDF; b) collisionless, so that the mean free path of the plasma particles is much longer than the largest characteristic 
scale length of the plasma; c) axisymmetric, so that the relevant dynamical variables characterizing the plasma (e.g., 
the fluid flelds) are independent of the azimuthal angle (p, when referred to a set of cylindrical coordinates {R, (p, z); 
d) acted on by both gravitational and EM fields. 

The kinetic formulation is intrinsically asymptotic. This means that the theory (in particular the GK theory 
formulated in the next section) is characterized by a suitable species-dependent dimensionless physical parameter 
eM,s = ^ 1) where s = i,e denotes the species index. Here rLs = v±ths/^cs is the species average Larmor 

1/2 

radius, with v±ths = {T±s/Ms} denoting the species thermal velocity perpendicular to the magnetic field and 
Q.CS = ZscB /MgC being the species Larmor frequency. Moreover, L is the characteristic length-scale of the spatial 
inhomogeneities of the EM field, defined as L Lb ~ Le, where Lb and Le arc the characteristic magnitudes 
of the gradients of the absolute values of the magnetic field B (x, t) and the electric field E (x, t) , defined as = 

max||^ln5 ,i = 1,3| and = max | -^InE ,i = 1,3|, where the vector x denotes x = {R,z). Then, in 

analogy with Paper I, we define a unique parameter em = max {sm.s, s = i, e}. For temperatures and magnetic flelds 
typical of AD plasmas, we have < Em <^ 1. 

In the following we will focus on solutions for the equilibrium magnetic field B which admit, at least locally, a 
family of nested and open axisymmetric toroidal magnetic surfaces {tp{ x)} = {ip{x) = const.}, where tp denotes 



Closed nested 
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surfaces 



Open nested 
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surfaces 



FIG. 1: Schematic comparison between the configuration of iocally closed magnetic surfaces considered in Paper I and the case 
of open magnetic surfaces analysed in the present study. 

the poloidal magnetic flux of B. See Fig. 1 for a schematic comparison between the configuration of locally closed 
magnetic surfaces considered in Paper I and the case of open magnetic surfaces analyzed in the present study. A set 
of magnetic coordinates ("0, 'J) can be defined locally, where ?? is a curvilinear angle-like coordinate on the magnetic 
surfaces ■(/'(x) = const. Each relevant physical quantity G(x, can then be conveniently expressed either in terms 
of the cylindrical coordinates or as a function of the magnetic coordinates, i.e. G(x, i) = G{ip,'0,t) , where the (f 
dependence has been suppressed due to the axisymmetry. 

We require the EM field to be slowly varying in time, i.e., to be of the form 

[E(x,4,i),B(x,4,i)], (1) 

with fc > 1 being a suitable integer. This time dependence is connected with either external sources or boundary 
conditions for the KDF. In particular, we shall assume that the magnetic field is of the form 

B = V X A = B"'^f{x,e'ljt) + B''=''{x,eijt), (2) 

where B*'^'-'' and B*^^* denote the self-generated magnetic field produced by the AD plasma and a finite external 
magnetic field produced by the central object (in the case of neutron stars or white dwarfs). For greater generality, 
we shall not prescribe any relative orderings between the various components of the total magnetic field, which are 
taken to be of the form 

B-'/ = /(x, eijt)V^ + V^p(x, eljt) x V^, (3) 
B^^* = V^AD(x,4,,t)x V(^. (4) 

In particular, here B^- = /(x, e\,ft)V(p and Bp = Vtpp{x, £\jt) x V(p are the toroidal and poloidal components of the 
self-field, while the external magnetic field B*^^* has to be purely poloidal, as a consequence of the axisymmetry, and 
is defined in terms of the vacuum potential '0D(x,e^jt). As a consequence, the magnetic field can also be written in 
the equivalent form 

B = /(x, e\it)V^ + V7^(x, e\,t) x V^, (5) 

where the function i/)(x, e^i) is defined as ^/'(x, e^i) = ^/'p(x, e^jt)-|-'(/'£)(x, e^i), with fc > 1 and (V', y^, i?) defining a set 
of local magnetic coordinates (as implied by the equation B • Vi/; = which is identically satisfied). Also, it is assumed 
that the charged particles of the plasma are subject to the action of effective EM potentials { 'I'^-'^''' (x, e'lft), A(x, e\ft)} , 
where A(x, e^.^i) is the vector potential corresponding to the magnetic field of Eq.®, while ^ {x, e'l.jt) is given by 

a>:^/(x,£^i) = a>(x,4,t) + |^i$g(x,4,o, (6) 

Zse 

with f {-x, e'lft) , $(x, e^^t) and <I>g(x, e^^) denoting the effective electrostatic potential and the electrostatic and 
generalized gravitational potentials (the latter, in principle, being produced both by the central object and the 
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accretion disk). The effective electric field E^-''-^ can then be defined as 



E-// = -V*'^-''^' - i— (7) 



III. GK THEORY FOR MAGNETIZED ACCRETION DISK PLASMAS 



In this section we recall the GK theory appropriate for the description of AD plasmas. Its formulation is in fact 
a prerequisite for the construction of the kinetic quasi-stationary equilibria to be developed later. The appropriate 
generalization of GK theory allowing for the presence of strong gravitational fields should in principle be based on a 
covariant formulation [see [23l426l |]. However, for non-relativistic plasmas within a gravitational field, the appropriate 
formulation can also be directly recovered via a suitable reformulation of the standard non-relativistic theory holding 
for magnetically confined plasmas [27l - l35| . 

In this case, the appropriate particle Lagrangian function can be represented in terms of the effective EM potentials 
|(j)e//(x, e^^t), A(x, e|ji)} , with fc > 1, where is defined in Eq.®. In terms of the hybrid variables z = (x,v) 
(with X and v denoting respectively the particle position and velocity vectors), this is expressed as 

/:,(z,^z,4,i) =r.P,-H3(z,4,i), (8) 

where = [M^v -I- %^A(x,e^^t)] and 

H.(z,4,0 - ^«2+Z,e$f /(x,4,i) (9) 

denotes the corresponding Hamiltonian function in hybrid variables. The GK treatment for the Lagrangian ([8]) 
involves the construction - in terms of a perturbative expansion determined by means of a power series in em - of a 
diffeomorphism of the form 

z= (r,v)^z' = (r',v'), (10) 

referred to as the GK transformation. Note that, in the following, we shall use a prime " ' " to denote a dynamical 
variable defined at the guiding- center position r' (or x' in axisymmetry) . Here, by definition, the transformed variables 
z' {GK state) are constructed so that their time derivatives to the relevant order in em have at least one ignorable 
coordinate (a suitably-defined gyrophase 0'). As an illustration, we show the formulation of the perturbative theory 
to leading-order in em- In this case the GK transformation becomes simply 

— r' _ w'xb' 

f^c, ' (11) 

where w' — w' cos (f)'e[ +w' sin ^'e'j , with (f)' denoting the gyrophase angle. In the following, the GK transformation will 
be performed on all phase-space variables z = (r, v), except for the azimuthal angle (p which is left unchanged and 
is therefore to be considered as one of the GK variables. Here b' — b(x', e^.^t), with b(x, e|^t)= B(x,e|^t)/i3(x, E^jt), 

while fl'^g — ^Ij^ and V^^j are respectively the guiding-center Larmor frequency and the effective drift velocity 

produced by E^'^-''-'^, namely 

V',fjiK,El,t)^^E';ff xh'. (12) 

The rest of the notation is standard, with u' and w' denoting respectively the parallel and perpendicular (guiding- 
center) velocities, both defined relative to the frame locally moving with velocity '^iff- follows that, when expressed 
in terms of the GK variables z', the GK Lagrangian and Hamiltonian functions, £(, and H'g, can be evaluated with 
the desired order of accuracy. In particular, to leading-order, i.e. neglecting corrections of 0(£^/) with n > 1, 

C'g — C's^^ + 0[em) and H'^ = H's^^ + 0{em), where £s and T-L's^^ recover the customary expressions 

4(i),r'.^A;-i^m;S'-H^W, (13) 
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^ m'^B' + ^ {u'h' + V^^^)' + Z,e$-, (14) 
with the magnetic moment m'^ ^ fi'^ = ^i^^, to leading order, while the gyrophase- independent modified EM potentials 



are 



C = <^'s", (15) 

M r 

a; = A' + -^(u'b' + V^,,), (16) 
Zjse 

in the same approximation. It is important to stress that the GK theory can be performed in principle to arbitrary 
order in em |27H35l | , thus permitting the explicit determination of m'^ and the modified EM potentials as well as the 
relevant guiding-center canonical momenta. 

A. First integrals of motion and guiding-center adiabatic invariants for AD plasmas 

The exact integrals of motion and the relevant adiabatic invariants corresponding respectively to Eqs.® and 
can be immediately recovered. By definition, an adiabatic invariant P of order n with respect to Em is conserved only 
in an asymptotic sense, i.e., in the sense that jp-gjlu-P = + 0{E^^f^), where n > is a suitable integer. First we 
notice that, under the assumptions of axisymmetry and of Eq.(IT]), the only first integral of motion is the canonical 
momentum p^pg = conjugate to the ignorable azimuthal angle ip: 

dip 

p^s ^ MsRv ■ Bu, + —^tp = ^^ij^s- (17) 
c c 

Since the azimuthal angle if is ignorable also for the GK Lagrangian C'g , it follows that the quantity p' = is an 

^ dip 

adiabatic invariant of the prescribed order, according to the accuracy of the GK transformation used to evaluate Cg. 
We shall refer to p'^^ as the guiding- center canonical momentum. In particular, correct to O(e^), with fc > 1, one 
obtains 



P^s = ^\ul + —V' , (18) 



which is an adiabatic invariant of 0(e^^^), with fc > 1. Furthermore, the total particle energy 

Es^^^v^+Zge<^lff{^,El,t), (19) 

with n > 1, and the GK Hamiltonian H'^ are also adiabatic invariants of order n. Finally, in GK theory, by con- 
struction, the momentum p'^,^ = dC'^/ dip' conjugate to the gyrophase, as well as the related magnetic moment m'^ 
defined as m'^ = jf^p'^ig, are adiabatic invariants. As shown by Kruskal (1962 [s^]) it is always possible to determine 
C'g so that m'g is an adiabatic invariant of arbitrary order in em- In particular, the leading-order approximation is 

™s = Ms = 2B' ■ 

Note that the allowance of slow time variations for Eg is an elementary consequence of assumption ([l} , which allows 
us to describe realistic configurations of AD plasmas which slowly evolve in time. 

B. Particle trapping phenomena 

GK theory permits explicit treatment of particle trapping corresponding to the existence of forbidden regions for the 
motion of charged particles arising from conservation of energy and magnetic moment. Conservation of the guiding- 
center Hamiltonian (|14p and the magnetic moment /i^ (leading-order approximation) give rise to some implications. 
Combining the two identities to express the parallel velocity u' , and using the definition ([T^. we find 




(20) 
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Therefore u' is a local function of the guiding-center position vector x' and, thanks to axisymmetry, of the corre- 
sponding flux coordinates {i})' Since the argument of the square root must be non-negative, this means that 
u' is only defined in the subset of the configuration space spanned by {4>',^') where this property holds. It follows 

that if the argument becomes null for given T-C's^'^'^ and ji'g, the parallel velocity must change sign so that the particle 
undergoes a spatial reflection. The points of the configuration space where this occurs are the so-called mirror points. 
The existence of these points may generate various kinetic phenomena in AD plasmas. In particular, for open mag- 
netic surfaces, particles can in principle experience zero, one or two reflections corresponding respectively to passing 
particles (PPs), bouncing particles (BPs) and trapped particles (TPs). In the present case, since the right hand side 
of Eg. ipO)) depends on the magnitude of the magnetic field (S'), the effective potential energy {ZsC^'f^^) and the 
centrifugal potential {^V^'jj), we shall refer to the TPs case as gravitational EM trapping. In Section VII we shall 
investigate some consequences of trapping phenomena for the dynamics of ADs. 



IV. CONSTRUCTION OF THE QSA-KDF: GENERALIZED SOLUTION 

In this section we show that the equilibrium generalized bi-Maxwellian solution for the KDF obtained in Paper I 
can be extended to QSA-KDFs describing axisymmetric AD plasmas with the following features: 

1) The KDF is also axisymmetric; 

2) Each species in the collisionless plasma is considered to be associated with a suitable set of sub-species (referring 
to the different populations mentioned above), each one having a different KDF; 

3) Temperature anisotropy: for all of the species, it is assumed that different parallel and perpendicular temperatures 
are allowed (with respect to the local direction of the magnetic field); 

4) Accretion flow velocity: a non- vanishing species dependent poloidal flow velocity is prescribed; 

5) Open, locally nested magnetic flux surfaces: the magnetic field is taken to allow quasi-stationary solutions with 
magnetic flux lines belonging to open and locally nested magnetic surfaces; 

6) Kinetic constraints: suitable functional dependencies are imposed so that the KDF is an adiabatic invariant; 

7) Analytic form: the solution is required to be asymptotically "close" to a local bi-Maxwellian in order topermit 
comparisons with previous literature dealing with Maxwellian or a bi-Maxwellian KDFs (see for example 0, H [SS] ) ■ 

Requirement 2) is suggested by observations of collisionless plasmas. For example, in the solar wind plasma both 
ion and electron species are described by superpositions of shifted bi-Maxwellian distributions. Requirements 1) - 7) 
clearly imply that the solution cannot generally be a Maxwellian. However, in analogy with Paper I, it is possible 
to show that they can be fulfilled b y a suitable modified bi-Maxwellian expressed solely in terms of first integrals of 
motion and adiabatic invariants (gI llOl [39.]. It follows that this is necessarily a QSA-KDF. A set of fluid equations 
can then readily be determined using this solution, expressed in terms of four moments of the KDF [corresponding to 
the species number density, flow velocity and the parallel and perpendicular temperatures]. These equations which, 
by construction, satisfy a kinetic closure condition, are also useful for comparing with previous fluid treatments. 

For consistency with the notation of Paper I, we again use the symbol "A " to denote physical quantities which 
refer to the treatment of anisotropic temperatures, unless otherwise specified, but in the present work, for greater 
generality, the symbol "*" is used to denote variables which depend on both the canonical momentum tjj^ts and the 
total particle energy Eg. 

In line with all of the previous requirements, it is possible to show that a particular solution for the QSA-KDF is 
given by: 

— /3« ^21) 



I I 

X exp <^ -— m^a^s 



which we refer to as the Generalized bi-Maxwellian KDF with parallel velocity perturbations. Here /,s is defined in 
the phase-space F = F,. x F^, where F^ and F^ are both identified with suitable subsets of the Euclidean space M.^ . 
The notation is as follows: 

ATs ^ 1^, (22) 



^ _ B' 

Q!*s = (23) 

At, 

K^s = Eg — £^s'!^*s, (24) 
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with Eg and ■0*s given by Eqs.([T9]) and (TTTI) respectively, while = J construction, i^^s has the 

dimensions of an angular momentum, while ro^^ has those of a frequency. In contrast with the solution obtained in 
Paper I, w^s is not necessarily associated here with a purely azimuthal leading-order velocity. In general Ks,s can, in 
fact, be represented as 

K^s ^ Es - p'^s^*s = - p'^sC^s- (25) 

Here H^s = Eg — ^^Tp,,s^*s has the same meaning as the analogous quantity used in Paper I, with 17*^ being related 
to the azimuthal rotational frequency. In Eq. ([25| is a frequency associated with the leading-order guiding-center 
canonical momentum p'^^ defined in Eq. (jl8p . which is an adiabatic invariant depending on u' and, by definition, is 
independent of the gyrophase angle. As we shall show at the end of this section, this feature can be used to require that 
the QSA-KDF carries a non-vanishing parallel flow velocity. This can be related to a net accretion flow arising in the 
AD plasma. Finally, by substituting Eg. pS)) into Eq. ((2T|) we reach the equivalent representation for the QSA-KDF: 



(2VM.)'/' {T\\.s) '' 



X exp 



T\\^, Til 



In order for the solution p6|) [or equivalently (|21l) ] to be a function of the integrals of motion and of the adiabatic 

invariants, the functions {A^^s} = |/3*s, Sli^, Tn^j, Sl^s, ^:,,s|, which we will refer to as structure functions, must be 

adiabatic invariants by themselves. To further generalize the solution of Paper I, we shall here retain a functional 
dependence on both the total particle energy and the canonical momentum, thus imposing the functional dependencies 

Ks{ip^s,Eg), (27) 

which will be referred to in the following as kinetic constraints. The kinetic constraints (j27p provide the most general 
solution for f^s- It can be shown that the physical motivation behind imposing these dependencies lies essentially in 
the fact that the asymptotic condition of small inverse aspect ratio (adopted previously in Paper I) is no longer valid. 
In the present context, the kinetic solution is no longer restricted to localized spatial domains in the disk but applies 
to the general configuration of open magnetic surfaces. This in turn implies that the structure functions are generally 
not simply flux- functions on the magnetic surfaces. In previous treatments (Paper I and Ref.js^), the structure 

functions were identified with 5^, Tn^s, r2*s| and |/3*s = -/V*s,T.s| for non-isotropic and isotropic generalized 

Maxwellian KDFs respectively. 
Some basic properties of f^^s are: 

Property 1: /*s is itself an adiabatic invariant, and is therefore an asymptotic solution of the stationary Vlasov 
equation, i.e., a QSA-KDF; 

Property 2: /^^ is only defined in the subset of phase-space where the adiabatic invariants p^^, H.'^^ and m'^ are 

defined. It follows that /^^ is suitable for describing both circulating and trapped particles (see the related discussion 
in Section 7); 

Property 3: all of the velocity-moment equations obtained from the Vlasov equation (and in particular the continuity 
and linear momentum fiuid equations) are identically satisfied in an asymptotic sense, i.e., neglecting corrections of 

Property 4: its velocity moments, to be identified with the fiuid fields, are unique once /^^ is prescribed in terms 
of the structure functions; 

Property 5: it generalizes the solution earlier presented in Paper I: a) by using both p'^^ and m'^ as adiabatic 
invariants and b) because of the new kinetic constraints. 

It follows immediately that the solution does indeed carry finite parallel velocity perturbations. Invoking the 
definitions ([T5)) and (P5|) . Eg. ([Ml) can be re- written as 



/3*sexp 



X exp ■ 
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where V,s = e^R^l^,s (V^*s,-£^s) and 



Here the function T' is defined as 



c 



(29) 



with J7|| 



(30) 



) 



,i;,). Note that U!, 



'\\*s ~ W^*'> \'r*s:-^sj- i""-""- I'iici.u is non-zero only if the toroidal magnetic field is non- vanishing. This 
quantity is independent of V*^ and is clearly associated with a parallel flow velocity (i.e., having both poloidal and 
toroidal components), referred to here as a parallel velocity perturbation. This perturbation enters the solution via 
the adiabatic invariant p'^^ and therefore its inclusion is consistent with the requirement that KDF is an adiabatic 
invariant. 

Finally we note that the same kinetic constraints ^7} also apply to the solution However, the functions 



/3*s exp 



f7|l^^ and T\\^s cannot be directly regarded as fluid fields, since they still depend on the single 



particle velocity via the canonical momentum t/j^.g and the particle energy Eg 



V. ANALYTICAL EXPANSION 

Based on Properties 1-5, in this section we determine an approximate analytical expression for f^,s obtained by means 
of suitable asymptotic expansions. These are carried out in terms of the following two dimensionless parameters: 



1) Ss- which is related to the canonical momentum ip^.s- This is defined as (cf Paper I) Ss 



where v^p = v • e. 



ip and Lips denotes the species particle angular momentum. We refer to the AD plasma as being 



W7 



i.e., it is the ratio between 



strongly magnetized if < 1; 

2) as', which is related to the total particle energy Eg. This is defined as as 

the kinetic energy and potential energy of the particle. For bound orbits Es < 0, and so cts < 1. 

In the following, we treat and as as infinitesimals of the same order, with Eg ^ as ^ and then Eg and ag can be 
used for performing a Taylor expansion of the implicit dependencies contained in the structure functions by setting 
^p^s = ip + O {Eg) and Eg = Zge<^'g^^ -t- O [ag) to leading order. This implies that the linear asymptotic expansion for 
the structure functions, obtained neglecting corrections of O [Egag) , as well as of O (e^) and O (cs ), with fc > 2, is 



A. 



A, + {i^^s - Tp) 



dJUs 



(Eg ~ ZsC^ 



•ff\ 



OA, 



dEg 



(31) 



where 



A., 



A 



(32) 



To perform the corresponding expansion for f^^g, we leave unchanged the dependence in terms of the guiding-center 
canonical momentum p'^g, while retaining the leading-order approximation for the magnetic moment only in the linear 
perturbation terms of Eg. pip . Then, it is straightforward to prove that for strongly magnetized and bound plasmas, 
the following relation holds to leading-order: 



f*s = Is {p'^g, m'g) [l + h\jg + h\,g] , 



(33) 



where h\)^ and lij^^ represent the so-called diamagnetic parts of /^^ (see the definition below). The definitions are 
then as follows: 
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First, the leading-order distribution fg {p'^st''^'s) expressed as 



X exp 



2 



A/.(v-V.-C/|^^b' 

ni „ - — 



2ri 



\\s 



(34) 



which we will here call the bi-Maxwellian KDF with parallel velocity perturbations. Here 
to the temperature anisotropy, the number density is defined as 



7^ 7?T- is related 

J is -Ills 



Us = rjs exp 



X, 



T, 



(35) 



and 



2 c 



(36) 



with Tjs denoting the pseudo- density. The function is defined as 



/2 



2a 



) 



(37) 



Note that = e^Rilg and J/^ = -gy^s define, respectively, the leading-order azimuthal flow velocity and the 
leading-order parallel velocity perturbation of the fluid. Then, the following kinetic constraints are implied from (|27p . 
to leading-order, for the structure functions: 



A, = A, (^, Z,e$f 



(38) 



Second, the diamagnetic parts h]^)^ and of /*s, due respectively to the expansions of the canonical momentum 
and the total energy, are given by 



r cMsR MsR,, I . ^ ^ 

{^[r, + y3] + ^r.^(v.e,). 



2Z,el Til, 



P'^sis I 2 

J li s 



(39) 
(40) 



Here Fi, i = 1, 5, is defined as 



Y3 
Y4 
Y5 



Y2 = flsil+Mss], 



Tn 



1 + ^C3. 



A2s 




Ik 


1 


T\\s 





tl'sCis 



(41) 
(42) 
(43) 

(44) 

(45) 



where Hg = Eg — ^^tps^s and the following definitions have been introduced: Ais = ^2s = —75;^, A^s 



d In T|| 



HnOa _ dal a _ In i r< — d\nf}s n — ^^n^lls n — dlnils — dSl r< — Sln^s 



0-0 



0*5 
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We should make a number of comments here: 

1) The functional forms of the leading-order number density, the parallel and azimuthal flow velocities and the 
temperatures carried by the bi-Maxwellian KDF, are naturally determined in terms of ijj and Zge^l^-^^ . The effective 
potential <&f-^^ is generally a function of the form ^^-^-^ — ^'^^ {k, e^jt) , with x = {R,z), since generally neither the 
gravitational potential nor the electrostatic potential are expected to be flux functions in the present case. Hence, in 
magnetic coordinates, it follows that the structure functions are of the form = (tp, {}, e|/i) ; 

2) The coefhcients Ais and Cis, i = 1,5, can be identified with effective thermodynamic forces: A^s carries the 
contribution of the parallel velocity perturbation, while the Ci^, * = 1,5, are due to the energy dependence contained 
in the structure functions; 

3) We stress that the energy dependence contained in the kinetic constraints is non trivial and cannot be included 
simply by redefining the structure functions (e.g., by transforming the magnetic coordinates). In fact, besides modi- 
fying the leading order structure functions (see point 1 above) , it gives rise to the new diamagnetic contribution . 
Ea.([33| is therefore a generalization of the analogous solution obtained in Paper I, which also appears in standard 
tokamak transport theory [s^, where the relevant structure functions were considered solely as flux functions. In- 
cluding the effect of the parallel velocity perturbations gives rise to contributions to hj^^ which are even with respect 
to u'; 

4) In the analytical expansion, we have assumed that the scale-length L is of the same order in as the characteristic 
scale-lengths associated with the structure functions; 

5) We have performed the analysis distinguishing between the different plasma species. Since this is an asymptotic 
estimation, the analytical expansion can be different for ions and electrons, particularly for the terms appearing in the 
diamagnetic part, depending on the relative magnitudes of the parameters Eg and a^- On the other hand, because of 
the double expansion and the energy dependence, the asymptotic solution for the two species can hold also in different 
spatial domains; 

6) The KDF (p'^g, also satisfies Property 2: namely, it is only defined in the subset of phase-space where the 
parallel velocity \u'\ is a real function. It is therefore suitable for properly describing particle trapping; 

7) Finally, we stress that the QSA-KDF (l26l) obtained here, reduces asymptotically to the expression reported in 
previous paper (see Eq.(lO) in Paper I) when the following conditions are satisfied: a) parallel velocity perturbations 
are ignored, namely the structure function is set to zero; b) closed nested magnetic surfaces are considered; c) 
large aspect ratio ordering, 1/(5 ^ 1, is invoked (see the definition in Paper I). In this case, the effective potential is 
solely a flux-function to leading order, while the diamagnetic contribution hjjg can be shown to be of higher order 
than /ik^. 



VI. MOMENT EQUATIONS 

In this section we discuss the connection between the kinetic treatment presented here and the corresponding fluid 
approach, obtained by describing the plasma in terms of a suitable set of fluid fields. The latter can in principle be 
specified as required by experimental observations and identified with the relevant physical observables. Important 
practical aspects of the present theory concern the explicit evaluation of the fluid fields associated with the QSA-KDF, 
and the conditions for validity of the relevant moment equations. 

For definiteness, let us require that: 

1. The KDF, the EM fields {E, B} and the corresponding EM potentials {$, A} are all exactly axisymmetric and, 
moreover, stationary in an asymptotic sense, i.e. neglecting corrections of 0(£^/^); 

2. The KDF is identified with the QSA-KDF f^g {Est4'*s, m's) which, by assumption, is required to be an adiabatic 
invariant of 0(e2/^)- By construction f„s (£'s, 'Jt-s) is a solution of the asymptotic Vlasov equation 

±-l^\nU = + O{e'l+'). (46) 

This equation holds by definition up to infinitesimals of O (^a/^), where n is an arbitrary positive integer; 

3. The magnetic field is taken to be of the form ([5]). 

As a basic consequence of these assumptions, the stationary fluid equations following from the Vlasov equation are 
necessarily all identically satisfied in an asymptotic sense, i.e., again neglecting corrections of O {s^m^)- In fact if 
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Z(x) is an arbitrary weight function, identified for example with Z = (l, v, u^) , then the generic moment of Eq. (|46)) 
is: 

^ d\Z^fi = + O{elt'), (47) 
where Tu denotes the appropriate velocity space of integration. Using the chain rule, this can be written as 

dip* dfts , dEs dfts , dm'^ df^s \ 



(48) 



/p^ ] dt dip^c ' dt dEs ' dt drn'^ 

= + {el+'). 

On the other hand, Eq. (|T7)) can also be represented as 



Z f*s 



d_ 

'di' 



Us-Z}=0 + O{sl+'), 



(49) 



which recovers the usual form of the velocity-moment equations in terms of suitable (and uniquely defined) fluid fields. 
For Z = (1, v) one obtains, in particular, that the species continuity and linear momentum fluid equations are satisfied 
identically up to infinitesimals of O (e^v/^) : 



(50) 



_^Vf xB = + O(e2/^). 



(51) 



Similarly, the law of conservation of the species total canonical momentum can be recovered by setting Z = f/^^s, 
namely 



d-^v — 
dt 



In the stationary case this implies the species angular momentum conservation law 



V 



for the species angular momentum 



(52) 



(53) 



(54) 



Here the notation is standard. In particular the following velocity moments of the QSA-KDF can be introduced: 

a) species number density 



d^vf^ 



b) species flow velocity 



c) species tensor pressure 



V* 



itot , 



d^vvf^ 



if* = ^ d\Ms (v - Vl"') (v - V*°*) 



(55) 



(56) 



(57) 
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d) species canonical toroidal momentum 

Llf = f cPv^^.Ss- (58) 
"'r„ c 

It is worth remarking here that the velocity moments are unique once the QSA-KDF [see Eq. (|2ip ] is prescribed in 
terms of the structure functions A*^,. On the other hand, as a result of Eqs. (l46l) and (j47p . it fonows that the stationary 
fluid moments calculated in terms of the QSA-KDF /^^ are identically solutions of the corresponding stationary fluid 
moment equations. In particular, imposing the quasi- neutrality condition in the sense 

Z^en':' = + (etf) (59) 
with A: > 2, the total fluid canonical toroidal momentum and the fluid angular momentum necessarily coincide, namely 

s—i.c s—i,e 

Let us now illustrate explicitly how it is possible to carry out such a calculation within the present theory. The 
evaluation of the previous fluid fields can be made by using the asymptotic analytical solution of the QSA-KDF ft^s 
derived in the previous section and given by Eq. p3p . For example, adopting this expansion in the limit of strongly 
magnetized plasmas, from Ea.(|55p the species number density becomes 

nf = ^ d^v[l[l + h],, + hl,]], (61) 

in which the diamagnetic corrections to the bi-Maxwellian KDF fs are polynomial functions of the particle velocity. 
Analogous expressions can also be obtained in a straightforward way for the remaining fluid moments. As pointed 
out in Paper I and subsequently in Ref. the expansion procedure for /,s can in principle be performed to higher 
order, allowing for the analytical computation of the corresponding quasi-stationary fluid fields and the determination 
of the relevant kinetic closure conditions for the stationary moment equations. In the present context we stress that 
the theory allows the treatment of multiple-species plasmas including, in particular, particle trapping phenomena. 
This is taken into account by proper definition of the velocity sub-space in which the integrations are performed. 
In fact, charged particles in both open and closed configurations can have mirror points (TPs and BPs) or be PPs, 
which are free to stream through the boundaries of the domain. These populations give different contributions to 
the relevant fluid fields and therefore require separate statistical treatments. The explicit calculation of fluid fields 
requires also a preliminary inverse transformation representing all quantities in terms of the actual particle positions 
(the FLR expansion, see Ea. ([TT|) '). This introduces further correction terms of order e\j, fc > 1, into the final 
analytical expressions. In contrast with the conclusion reached in Paper I, here we expect these FLR corrections to 
be non-negligible due to the requirement £m,s ^ £« holding for open-field configurations. 



VII. SLOW TIME-EVOLUTION OF THE AXISYMMETRIC QSA-KDF 

In this section we investigate the temporal evolution of the axisymmetric QSA-KDF, consistent with the assumptions 
of Section 2 and the results of Sections 3 and 4. Two different issues must be addressed: giving an estimate of 
the maximum time interval over which the QSA-KDF can be regarded as an asymptotic stationary solution; and 
determining the solution of the Vlasov equation for time intervals longer than the equilibrium one. 

For our explicit determination of the time evolution of the QSA-KDF, we make the following assumptions: 

1) That the plasma can be treated as a continuous medium in the kinetic description. This requires that the 
species kinetic equation holds on time and spatial scales which are much longer than the corresponding Langmuir 
characteristic times and Debye lengths; 

2) That we are considering timescales much shorter than the species characteristic collisional time tc so that it is 
appropriate to use the Vlasov equation; 

3) That the species KDF and the EM fields vary slowly in time and space with respect to the corresponding Larmor 
times and radii, so that the GK description is valid; 
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4) That the EM and gravitational fields vary slowly in time, in accordance with Eq.(IT]), so that the total energy Eg 
is an adiabatic invariant. In particular, we require: 

which implies that Tis^lni^s ~ O ^^jf/j^ j , with n > 0. Consistently with the properties of solution (pS)) . we take 
n = as a specific case. Note that from here on, tls = ^p- will denote the species characteristic time associated with 
the Larmor rotation (the Larmor rotation time). Since ^f^*s — it follows that 



dt"*s dt dEs 

C Otj. 



(63) 



5) That the magnetic moment m'^ and the guiding-center canonical momentum p'^^ can be taken as adiabatic 
invariants of O (^e\j ^ , with j > n. The ordering tl^-^ lnp!^s ^ O (e^,^ ^) holds for the leading-order expression for 
p'^g adopted here as follows from Eq. p^ and the fact that, by definition, higher-order correction terms, Ap'^^, to p'^^ 



are independent of the gyrophase angle </>'. In fact, denoting by C's"^' the second-order GK Lagrangian, Ap' can be 



estimated as Ap' = — L's"^^ — C's'^^ where, by construction, and C's'^' are both gyrophase independent. Note 

^ dip L J 

that the assumption made here requires the construction of a higher-order GK theory in order to correctly determine 
m'j, to the required order in the Larmor-radius expansion. 

The time evolution of the QSA-KDF is in principle determined by two different mechanisms: the explicit time 
variation of the EM and gravitational fields, and the time variation of the guiding-center adiabatic invariants. However, 
the choice of the orderings in 4) and 5) above, allows the time dependence produced only by the EM and gravitational 
fields to be singled out. 

When assumptions 1) - 5) above hold, it follows that t^s^ ln/*s = + , with ti > being determined by 

Ea. (|62p . Then, ignoring higher-order corrections 

^In/T.-^^s, (64) 



'(1) ^^A ^'(2) 



where 



dt - dt 



dEs ' dEs Ty,, dEs 

H,s 1 , p'ipsC*s\ 9 In Til, 



T\\^:s 2 T\\^s J dEs 

1 / Zse dn,s\ 

^-U^ ' (65) 



T\\^s V c dEs , 

and so the solution /h.^ can be regarded as an exact kinetic equilibrium for all times t > such that 

TLs < t < tsnp < Tc, (66) 

where tgup = -^ft- Within the scope of the above assumptions, we now determine the dynamical evolution equation 
which describes the slow time-evolution of the QSA-KDF /*s, for time intervals such that t is within 

tsup <i<Tc. (67) 

In analogy with Ref.fsHi, we denote by 

fs-Us + g's (68) 

the exact solution of the collisionless Vlasov equation, for which = 0. Here g'g is referred to as the reduced 
KDF. Following the discussion in Ref.js^, regarding the evaluation of ^g^: it is straightforward to prove that g's is 
gyrophase independent, to lowest order, in the sense that ^% = 0. Therefore, identifying the GK variables with the set 
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FIG. 2: Schematic view of the configuration geometry (not to scale) and meaning of the notation. 



z = ^1?', ip,p'^g,'H's^\m'g, 0'^ , we shall assume that g'^ is axisymmetric and of the form g'^ = g'^ {d' ,p'^^,'H's^\m'g,t^ . 

The gyro-averaged dynamical equation for g'^ can then be obtained to next order by introducing the gyro-average 
operator {■■■)^, defined as 

(69) 

with the operation being performed while all of the other GK variables are held fixed [s^. It follows that, to leading- 
order, -^g'g = -^g'^ + d ^7(7^, where the time variation of the guiding-center magnetic coordinate d' is given by 



■ V'l?' to leading-order, with the equation of motion for r following from the gyrokinetic 

Lagrangian [e.g. from the leading-order Eq. (jl3|) ]. Then, consistently with these assumptions and ignoring higher-order 
corrections, it is found that the GK reduced KDF g'^ obeys the reduced GK- Vlasov equation 

d-t's+^g^9s = -{f*sSs^)^^. (70) 



It follows that, to leading-order 



Denoting/^ = [^^s,H^s-,v'^s-,'^'s) ■• f's is then defined as = Fs (^-^p'^s^'^'s^\p'^s^'m'^- The remaining gyrophase 
average in the last equation can be performed in a straightforward way using Eas. ([62)) and (|65ll . 

Eg. ([70)1 clearly also holds in the time interval ([M)) . and so it determines the slow time-evolution for all times 
tls t <ti tc- For consistency, the non-stationary Maxwell equations must also be solved with the same accuracy. 
Eg. (1701) must be supplemented by appropriate boundary conditions: for open magnetic surfaces with boundaries 
prescribed on a given magnetic surface ip — const., at 1) — -di and 1) — 'd2, with i?i < 1)2 and z^i, '&2 representing the 
internal and external boundaries, these are defined respectively either by prescribing fs (i?i) = /i^^ or fs (1^2) = /i^"* 
(see Fig. 2 for a schematic view of the configuration geometry and the meaning of the notation). Both /i^-* and fs"^^ 



are necessarily of the form (|68p but their moments remain arbitrary in principle. As indicated below, this is essential 
for making comparisons with experimental observations. 

The results obtained here have important consequences for the kinetic description of slow time-evolution of colli- 
sionless AD plasmas. Fluid fields and moment equations can be explicitly determined in terms of Eg. (|70|) by invoking 
the perturbative expansion outlined in Section 5 and the relations given in Section 6. 
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VIII. THE AMPERE EQUATION AND THE KINETIC DYNAMO 



In this section we apply the kinetic solution for the QSA-KDF to discuss the properties of the Ampere equation 
and the implications for the self-generation of magnetic field by the quasi-stationary AD collisionless plasma. We 
refer here to this phenomenon as a quasi- stationary kinetic dynamo effect. Generalizing the treatment presented in 
Paper I, the Ampere equation for the self magnetic field becomes: 

VxB-'^ = — (J^+J^+J^), (72) 

where B*'^'-'' has been defined in Eq.Q and here we have distinguished between the contributions arising from PPs, 
BPs and TPs, denoting the corresponding total current densities as J"^, and J^. As described in Sections 5 and 6, 
these fluid fields can be calculated in closed analytic form to the required order, by using the asymptotic expansion 
of the QSA-KDF. This gives: 

s—i,e s—i,e 

- Y Zse I d?vw [I [1 + h],, + /i^,] } (73) 

s—i,e 

for I ~ T, B, P and where Pj^ denotes the appropriate velocity space domain of integration for trapped, bouncing and 
passing particles respectively. For convenience of notation, in the following we shall denote as J = J"^ + -I- J''^ the 
total current density entering Eg. ([7^. It is possible to prove that in the case of open magnetic surfaces the total 
current density J in general has non-vanishing components along all of the three directions identified by the set of 
magnetic coordinates ("0, (/?, i?). Hence, J can be represented as 

J = ( J^.Vi? X Vip, J^Vip, J^iVtA X Vip) . (74) 

Let us now proceed with the study of the Ampere equation. The toroidal component of Eg. ([7^ gives, as usual, the 
generalized Grad-Shafranov equation for the poloidal flux function ipp: 

47r 

AV^ = - — (75) 

where the elliptic operator A* is defined as A* = i?^V • (i?~^V) . The remaining terms of Eq. ([72l) along the directions 
V'i? X V(p and x Vip give two equations for the toroidal component of the magnetic field I/R. These are respectively 

'^J., (76) 



dil) c 
dl An 



yielding the constraint 



J^, (77) 
c 



(78) 



which is a solubility condition for the structure functions. In this regard we notice that as a consequence of the kinetic 
constraints the function / in the previous equations is of the form I {ip , "ff , e'l^t) , i.e., in contrast to Paper I it is no 
longer a flux-function. Therefore, the solubility condition ([78]) can always be satisfled. Eqs. (|75| -([78 |) therefore provide 
consistent solutions for both poloidal and toroidal self magnetic flelds in a collisionless AD plasma. 

It is remarkable that in principle all of the populations of charged particles (PPs, BPs and TPs) can contribute to 
the generation of the toroidal magnetic field. More precisely, the following mechanisms can be involved: 

^1) FLR and diamagnetic effects, driven by temperature anisotropy, of the same kind as those described in Paper 

I; 

#2) Parallel velocity perturbations U'^^^, which generate a poloidal flow velocity, giving a related contribution to 
the electric current density through and J^; 

#3) FLR effects driven by the remaining thermodynamic forces (see Section 5). These contributions are produced 
by the diamagnetic KDF and arise because of the asymptotic ordering introduced here; 
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#:4) Gyrophase-dependent contributions driven by the same thermodynamic forces. These are originated by the 
inverse GK transformation of the guiding-center quantities in the QSA-KDF. 

As discussed above, contributions #2 and #4 were neghgible under the circumstances discussed in Paper I. Therefore 
they should be considered as characteristic features of open-field configurations. 

We refer to the mechanism of self-generation of both poloidal and toroidal magnetic fields as a quasi- stationary 
kinetic dynamo effect. In contrast to customary MHD treatments, this type of dynamo effect occurs in the absence of 
possible instabilities or turbulence phenomena. In particular, in the case of TPs, the self generation of toroidal field 
could take place even without any net accretion in the domain of interest, in presence of open magnetic field lines. 
This phenomenon is analogous to that treated in Paper I for closed-field configurations. In particular, the toroidal 
field is associated with the existence of torques which cause redistribution of angular momentum, producing radial 
infiows and outfiows of disk material. As a consequence, various scenarios can be envisaged in which stationary radial 
fiows and kinetic dynamos are present in AD plasmas, both affected by processes of type ^1-^4. 



IX. QUASI-STATIONARY ACCRETION FLOW 

Let us now consider specifically the application of the kinetic solution developed here to the investigation of the 
accretion process in AD plasmas. 

The inward accretion fiow in ADs is usually "slow" in comparison with the characteristic Larmor time tls- For 
example, AD plasmas with B ^ 10^ — lO^G have Hydrogen-ion Larmor rotation times in the range r^i ~ 10"'' — 10~^^s 
which is shorter than the dynamical timescale at most relevant radii. For typical plasma densities and temperatures 
in the range rii ~ 10^ — lO^^cm'^ and Ti ^ 1 — lOkeV, the (Spitzer) ion collision time (below which the plasma can 
be considered coUisionless) is in the range tc ~ 10^ — lO^s (the upper value corresponding to high temperature and 
low density). Independent of the physical origin of the accretion process, we can therefore expect that the present 
theory correctly describes phenomena occurring on all time-scales in the range r^i < t < tq. 

We next determine the local poloidal and radial fiow velocities for the various particle sub-species. By definition, 
these are given by 

^ps — ^ s ' 

= E i / d'vW-e,]f:,[l+gi], (79) 

sub— species u 

Vns^Vs-en= ^ -^4,, (80) 

sub— species 

4, EE / d\[veR]U[l+gi], (81) 

where Bp = ^Vix Vy| ^^'^ ^ ^^'^ summations are performed over the particle sub-species for I ^ T, B, P. 
We stress that the velocity-space integrals indicated above must contain the contributions from PPs, BPs and TPs 
and so Jrs = JJis + Jj^s + -^fls' where J]^^, J^^ and J^^ are the corresponding mass currents. As an example, let us 
consider the leading-order contributions obtained ignoring FLR corrections. Explicit calculation gives 

Vps = C/pb-Gp, (82) 
Vr, ^ [/pb-Gfl, (83) 

where U\\s = U\\s (V')''^)£mO = functional dependence = {i^ t t is prescribed by the kinetic 

constraints (p8| . We stress that the precise form of (il^ , 1} , e'ljt) still has to be chosen to satisfy the solubility 
constraints imposed by Ampere's law (see the discussion in previous section) and so the radial mass current density 
Jrs is generally a function of the form Jrs = Jrs (x, £4/*) — dRs{'4',^,£M^)- We are interested in situations where 
there is a net radial accretion flow i.e. where the average radial mass current {{Jrs)) = z2-zi Izi '^Rs'^^ (with zi and 
Z2 being suitably prescribed) is negative. There are local contributions to {(Jrs)) from TPs, BPs and PPs, but the 
overall accretion flow is mainly associated with PPs. 

Let us show that such a solution exists. We seek particular Vlasov-Maxwell equilibria which are globally quasi- 
neutral, in the sense of Eg. (15^1) . These equilibria are uniquely defined once $(x, e^t), tp (x, e^t), I {tp , 1} , e'l^t) and 
the structure functions are prescribed. The latter, by definition, are arbitrary smooth real functions of the specified 
variables as required by the kinetic constraints. Notice that, if the quasi-neutrality condition is valid, an analytical 
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solution for the ES potential can be obtained as shown in Paper I. Furthermore, we consider an example in which by 
assumption 

1. C/||s is non- vanishing and dominant with respect to FLR effects, so that Eqs. (l82|) and (l83l) apply. 

2. Particular solutions have a definite parity property with respect to the spatial reflection z ^ —z. As a specific 
case, the poloidal fiux ^ is assumed here to be antisymmetric, i.e., ip (z) = —ip (— z), while both the toroidal and 
poloidal magnetic fields are symmetric. As a consequence, the toroidal current density must be antisymmetric. 
This can be realized only if a vertical electric field is present (i.e., one in the z direction), consistent with the 
quasi-neutrality condition. 

If these assumptions are valid. Ampere's law demands that, to leading order, i.e., neglecting diamagnetic FLR 
effects, 

- ^^Y.ZsenTU\\s. (84) 



dij} c 
dl 



d-d 



0, (85) 



namely 1 = 1 {ij)) at this order of approximation. Therefore, in this case a solution consistent with the requirement 
of net radial accretion flow and Vlasov-Maxwell equilibrium is obtained imposing that the species number density 
ri*°* is even in ?/;, while the species structure function is odd with respect to the same variable. A solution of this 
type is consistent with the angular momentum conservation law (1531) ; in order to obtain the solution, suitable kinetic 
boundary conditions must be prescribed (see the discussion following Eq. (l7ip ). This proves that stationary accretion 

solutions exist and are admitted by the present kinetic theory for the "incoming" QSA-KDF, namely for f^g in the 

subset V ■ < (see Fig. 2). The same conclusion is in principle applicable for outflows, by appropriate prescription 

of the "outgoing" QSA-KDF /,s in the subset v • > 0. In fact, the angular momentum conservation law ([55)1 

allows both inward and outward radial fluid velocities for each species, namely having V*°* -efl < or > respectively. 
Indeed, for a collisionless plasma the species tensor pressure is generally non isotropic (see the related discussions 
in Paper I and Ref . [l6j ) such that Eq. ((53| is identically satisfied. Unlike the customary view based on ideal MHD, 
for which a self-consistent treatment of inflow and outflow solutions is usually difficult, within the present theory 
both infiows and outfiows can occur independently and are described consistently by their respective QSA-KDFs. In 
particular. Eg. ((55)) shows that radial flows arise due both to the parallel velocities U\\s and to the kinetic effects carried 
by the FLR diamagnetic corrections. As a result, species radial fiow velocities appear necessarily in combination with 
non-isotropic tensor pressures and a non-vanishing toroidal magnetic field. In conclusion, the theory predicts the 
possibility of having purely infiowing matter in quasi-stationary AD plasmas, or of having co-existing inflows and 
outflows. 



X. CONCLUSIONS 



In this paper, a consistent theoretical investigation of the slow kinetic dynamics of collisionless non-relativistic and 
axisymmetric AD plasmas has been presented. The formulation is based on a kinetic approach developed within 
the framework of the Vlasov-Maxwell description. We have considered here plasmas immersed in quasi-stationary 
magnetic flelds characterized by open nested magnetic surfaces. This can be appropriate for radiatively inefficient 
accretion flows onto black holes, some of which are believed to be associated with a plasma of collisionless ions and 
electrons having different temperatures, and there can be other related applications to the inner regions of accretion 
flows onto magnetized neutron stars and white dwarfs. The discussion presented here provides a background for future 
investigations of instabilities and turbulence occurring in these plasmas. 

We have shown that a new type of asymptotic kinetic equilibria exists, which can be described by QSA-KDFs 
expressed in terms of generalized bi-Maxwellian distributions. These solutions permit the consistent treatment of a 
number of physical properties characteristic of collisionless plasmas. The existence of these equilibrium solutions has 
been shown to be warranted by imposing suitable kinetic constraints for the structure functions entering the definition 
of the QSA-KDFs. In terms of these solutions, the slow dynamics of collisionless AD plasmas has been described 
by means of a suitable reduced GK-Vlasov equation. In addition, the theory permits the consistent treatment of 
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gravitational EM particle trapping phenomena, allowing one to distinguish between different populations of charged 
particles. 

We have shown that the kinetic approach is suitable for the description of quasi-stationary AD plasmas subject 
to accretion flows and kinetic dynamo effects responsible for the self-generation of both poloidal and toroidal mag- 
netic fields. Four intrinsically-kinetic physical mechanisms have been included in the treatment of this, related to 
temperature anisotropy, parallel velocity perturbations and FLR-diamagnetic effects. 

The novelty of the present approach, with respect to traditional fluid treatments, lies in the possibility of explicitly 
constructing asymptotic solutions for the fluid equations: the calculation of all of the relevant fluid fields involved 
(e.g. the plasma charge and mass current densities and the radial flow velocity) can be performed in a straightforward 
way using a species-dependent asymptotic expansion of the QSA-KDF. 

We believe that this study makes a relevant contribution for the description of two-temperature collisionless AD 
plasmas and the improvement of our understanding of their physical properties. The kinetic treatment developed 
here can also provide a convenient starting point for making a kinetic stability analysis of these plasmas. 
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